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Abstract 

We study the interband Landau-Zener tunneling of a quantum particle in the Hall configuration, 
i.e., in the presence of normal to the lattice plane gauge field (for example, magnetic field for a 
charged particle) and in-plane potential field (electric field for a charged particle). The interband 
tunneling is induced by the potential field and for vanishing gauge field is described by the common 
Landau-Zener theory. We generalize this theory for non-zero gauge field. The depletion rates of 
low-energy bands are calculated by using semi-analytical method of the truncated Floquet matrix. 
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I. INTRODUCTION 



In recent years we have seen a recovering of interest to the Landau-Zener tunneling 
(LZ-tunneling) in periodic structures. Although this phenomenon was originally discussed 
with respect to Bloch oscillations (BO) of crystal electrons in a strong electric field 1-4], 
nowadays the most successful experimental systems are semiconductor superlattices 5ruj], 
one-demensional arrays of optical waveguides pHl 1| . and cold atoms in quasi ID optical lat- 
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ow access to many different aspects of LZ-tunneling including 
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26], LZ-tunneling in binary (double-periodic) lat- 



25], time- resolved LZ-tunneling 



nonlinearity caused by atom-atom interactions 
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26|, etc. 



The above cited papers refer to effectively one-dimensional systems. Another direction 



of research is BO and LZ-tunneling in two-dimensional periodic structures 
parison with ID lattices LZ-tunelling in 2D latices depends not only on the field magnitude 
but also on the lattice geometry (square, hexagonal, etc.), direction of the field vector with 
respect to the lattice primary axes, and particular properties (for example, separability) 
of the periodic potential. In the present work we address LZ-tunneling in 2D lattices in 
the presence of normal to the lattice plane magnetic field (electron systems) or artificial 
gauge field (cold atoms and twisted waveguide arrays), which mimic the magnetic field for 



charge neutral particles 33j, I35H371]. Using the solid state terminology we shall refer to 



these systems as quantum particle in the Hall configuration. In what follows we define the 
notion of LZ-tunneling for the quantum particle in the Hall configuration and obtain an 
estimate on its rate. This estimate is highly demanded to formulate the validity condition 
of the tight-binding approximation, that is widely used in physical literature to analyze the 
cyclotron-Bloch dynamics of the quantum particle in the Hall configuration 
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Mi- 



ll. THE SYSTEM 



Let us consider a square lattice, where the electric field is aligned with the y axis. Then, 
using the Landau gauge for a magnetic field the dimensionless Hamiltonian of the quantum 
particle in the Hall configuration reads, 

1 



H 



p 2 x + (p y -Bxy +V(x,y) + Fy 



(1) 
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FIG. 1: (color online). Survival probability as the function of time for F = 0.05 and B = (dashed 
line) and B = 1/16-7T ~ 0.02 (solid line). The lattice parameters are v x = 0.5 (J x = 0.0431) and 
v y = 0.25 (J y = 0.0741). The time is measured in units of the Bloch period Tb = 1/F. 

where B and F are the magnitudes of magnetic and electric fields, respectively. For vanish- 
ing electric and magnetic fields the spectrum of the Hamiltonian (CQ) consists of the ground 
Bloch band separated from the rest of the spectrum by a finite energy gap. We assume 
that the initial state of the particle belongs to this ground band. A finite F induces Bloch 
oscillations in the band and simultaneously causes LZ-tunneling across the energy gap. In 
what follows we shall consider the separable periodic potential, 

V(x, y) = v x cos x + v y cos y 

(the lattice period is set to 2tt). In this case the rate of tunneling can be readily calculated 
because the 2D Hamiltonian factorizes into two ID Hamiltonians if B — 0: 

H x = y +v x cosx , (2) 
— p 2 

Hy = -J + v y cosy + Fy . (3) 

Thus we can use the known results for ID lattices - population of the ground Bloch band de- 
creases exponentially in time with the rate T given by the inverse lifetime (resonance width) 



of the ground Wannier-Stark states [161 ] . Ignoring the phenomenon of resonant tunneling 
the F dependence of the rate Y is given by the celebrated Landau-Zener equation, 

f(F)~Fexpf-^) , (4) 

where h is proportional to the square of the energy gap A y of the Hamiltonian For 
the purpose of future comparison the dashed line in Fig. [1] shows population of the ground 
Bloch band as the function of time for v x = 0.5, v y = 0.25, and F = 0.05. The initial 
wave-function corresponds to the ground state of the Hamiltonian (JT|) for vanishing electric 
and magnetic fields, i.e., to the Bloch wave with zero quasimomentum. The steps in P(t) 
occur when the quasimomentum n y , which evolves according to the linear law K y = Ft, 
crosses the boundary of the Brillouin zone. 



III. MAGNETIC BANDS 



Before proceeding to LZ-tunneling for B ^ we need to discuss the notion of ground 
magnetic bands. By those we mean magnetic bands that are originated from the ground 
Bloch band. One easily finds them by using the tight-binding approximation to the original 
Hamiltonian, 

H tb = E J2 \l, m) (I, m\-^ £ (|Z + 1, m) (I, m\ + h.c.)-^- £ (|Z, m + 1) (/, m|e M + h.c.) , 

l,m l,m l,m 

(5) 

where \l, m) denote the Wannier functions associated with the ground Bloch band, is the 
on-site energy, J x>y are the hopping matrix elements, and a = 2itB is the Peierls phase. If 
a = the spectrum of ([5]) is given by 

E(k x , K y ) = E — J x cos k x — Jy COS Ky 

and the eigenfunctions are Bloch waves, ~ Hi, m exp[z(K x Z + K y m)\l,m). If a ^ so- 
lutions of the stationary Schrodinger equation with the Hamiltonian have the form 
l^) ~ J2i,m exp(zK J/ m)6/|Z, m), where the coefficients b[ satisfy the Aubry-Andre equation, 

- + h-x) - J y cos(27raZ + K y )bi = Eb x , (6) 

which for J x = J y coincides with the Harper equation. If a is a rational number we can apply 
the Bloch theorem and, hence, eigenfunctions of (j6]) are labeled by the quasimomentum k x 
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FIG. 2: (color online). Ground Bloch band for (J x , J y ) = (0.0431, 0.0741), top, and magnetic bands 
originated from this band if a = 1/8, bottom. 

defined in the reduced Brillouin zone. Thus for rational a = r/q the ground Bloch band splits 
into q magnetic bands. As an example the lower panel in Fig. [2] shows the ground magnetic 
bands for a = 1/8 and (J x , J y ) = (0.0431, 0.0741), which corresponds to (v x , v y ) = (0.5, 0.25) 
in the continuous Hamiltonian. Notice that in Fig. [2] we use the extended Brillouin zone 
picture. We also mention that in the considered case J y > J x the magnetic bands are flat 
in the k t direction. 
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For deep lattices the above tight-binding approach provides a reasonable approximation 
to the ground magnetic bands. However, to address LZ-tunneling we need to find not only 
ground bands but also the spectrum above the energy gap. In the rest of this section we 
discuss a different method of calculating the ground magnetic bands which accomplishes 
this task. 

Let us assume that v x is large enough to justifies the tight-binding approximation in the 
x direction. Then we can use the ansatz 



l=— oo 



(7) 



where 4>i(x) are the ID Wannier functions associated with the ground Bloch band of the 
Hamiltonian (j2J) . Substituting (JTJ) into the stationary Schrodinger equation with the Hamil- 
tonian (JTJ), where we temporally set F = 0, we have 



where 



-(py-al) +v y cosy 



(8) 



(9) 



The eigenf unctions of the Hamiltonians (Q are Bloch waves with shifted dispersion relation. 
Namely, if E(k) is the Bloch spectrum of the Hamiltonian H^°\ then for I ^ we have 



E {1 \k) = E^in + al) , 
see Fig. [3]^a). Next, using the Fourier expansion for the Bloch wave, 



^ {l \y) = exp(wy) J2 c® (k) exp(iny) 



(Oi 



(10) 



(11) 



the system of partial differential equations flS} reduces to the system of algebraic equations 
for the coefficients cjp: 



J X 

T 



\{n + K + alfc^ + V -l 



c (0 + c {0 



ES . 



(12) 



In general case of arbitrary a the index I in (151 fT2"j) runs from minus to plus infinity. However, 
if a = r/q is a rational number we can restrict I to one magnetic period, 1 < I < q. In this 
case Eq. ffl2|) should be accomplished by the periodic boundary conditions 



-,(9+1) 



,(1) 



(13) 
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FIG. 3: Magnetic bands for a = 1/8 and (J x ,v y ) 
(0.0431,0.25), right panel. 



(0.0,0.25), left panel, and (J x ,v y ) 



The system of algebraic equations (|T2|) together with the boundary condition ( TT3|) provide 
an alternative method for calculating the ground magnetic bands. 

The right panel in Fig. [3] shows the solution of Eqs. (I12|ll3p for a — 1/8, v y — 0.25, and 
J x = 0.0431, which corresponds to v x = 0.5. Magnetic bands are plotted as functions of 
the quasimomentum k = K y for a single value of the quasimomentum k x = 45|. This 
figure should be compared with Fig. Mjo) showing the magnetic bands in the tight-binding 
approximation. We note that in Fig. [3] we intentionally restricted the upper limit of the 
energy axis to a relatively low value because for higher energies the second Bloch band 
of the Hamiltonian H x contributes the spectrum. However, for our aim of studying LZ- 
tunneling it is sufficient to have a fragment of the actual energy spectrum just above the 
energy gap. 



IV. LZ-TUNNELING 



The structure of the eigenvalue equation provides insight into physics of LZ-tunnelinj 
in the presence of a magnetic field. To address this phenomenon the ID Hamiltonians ([9j 
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should be accomplished by the term Fy and we should solve the non-stationary Schrodinger 
equation instead of the stationary one. Thus the original 2D problem reduces to q coupled 
ID Landau-Zenner problems. Below we analyze the effect of this coupling on the tunneling 
dynamics by using two different approaches. 

In our first approach we numerically solve the time-dependent Schrodinger equation and 
project the solution on the subspace of ground magnetic bands. In more detail, we solve the 
following differential equations for the coefficients c®, 



where the quasimomentum k' linearly depends on time as k' = K+Ft. [This equation follows 
from (11 ip where one substitutes k' instead of k.\ The initial wave-fuction corresponds to one 
of multi-degenerate ground states of the system for F — 0. Thus initially only the lowest 
magnetic band is populated. Depending on the electric field magnitude F, we observed 
three different regimes of the tunneling dynamics. For very small F population of the lowest 
magnetic band stays close to unity during the whole simulation time. If F exceeds some 
critical value F cr , all q ground magnetic bands become involved into dynamics yet their 



+ \{n + k' + alfcU + V f (H + cti) , (14) 



populations sum up to unity with high accuracy |46|. Finally, for even larger F we observe a 
decrease in the total population of ground magnetic bands, see Fig. [TJ Comparing P(t) for 
a = 1/8 (solid line) and a = (dashed line) it is seen that finite magnetic field decreases the 
rate of LZ-tunneling and smoothes oscillations of P(t), although the overall decay remains 
exponential [see inset in Fig. [TJ. 

In the above simulations the initial wave function belongs to the lowest magnetic band. In 
fact, the rate of tunneling across the energy gap strongly depends on which magnetic band is 



initially populated. To find decay rates T ~ o 



the individual magnetic bands we employ the 



truncated Floquet matrix method of Ref . [28] , adopted to the currently considered problem. 



Namely, using the ansatz (iTTj) and the Schrodinger equation (fl4j) we calculate the (formally 
infinite) matrix of the evolution operator over one Bloch period Tg = 1/F and truncate it to 
a finite size. Notice that when calculating the Floquet matrix we explicitly use the periodic 
boundary conditions ( fl3l) . Thus the matrix is truncated only with respect to the index n, 
n < N . (The method rapidly converges if N is increased, in our calculations we use N — 7.) 
The eigenvalues Xj of the truncated Floquet matrix are known to be the complex poles of 
the scattering matrix. Then the individual decay rates Tj of ground magnetic bands are 
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FIG. 4: Decay rate T of the ground Bloch band (B = 0), left panel, and decay rates Tj of the 
ground magnetic bands for a = 1/8 (B = a/2n), right panel, as functions of the inverse electric 
field magnitude. 

found from the equation |Aj| 2 = exp(— T^Tb), where Xj are the first q eigenvalues which are 
closest to the unit circle. 

The decay rates Tj for a = r/q = 1/8 are shown in Fig. H^b) as functions of the inverse 
electric field magnitude. The right panel of this figure should be compared with the left 
panel, which shows the decay rate of the ground Bloch band for B = 0. As mentioned in 
Sec. mi for B = the functional dependence of the decay rate is approximately given by 
the Landau-Zener formula (j4j), where deviations are due to phenomenon of the resonant 
tunneling 16[. It is seen in Fig. |4]that the resonant tunneling also takes place if B ^ 0. 
Thus we can decompose Tj into two terms, 



r,(F) = f,(F) + rf(F) 



(15) 



where fj(F) ~ Fexp(—bj/F) and T^(F) is the oscillating part. Analyzing the values of the 
coefficients bj we conclude that LZ-tunneling is suppressed for lower magnetic bands, j q, 
but enhanced for higher bands, j ~ q. Notice that this effect is well pronounced only for 
weak electric fields, while in the strong field regime the decay rates Tj are essentially the 
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same and approximately coincide with that for B = 0. 

To conclude this section we briefly discuss the case of irrational a. Since in the wave- 
function simulations we do not use any boundary conditions, we can directly compare the 
survival probability P(t) for rational and irrational a. However, this approach gives reliable 
results only for strong electric fields, where P(t) essentially differs from unity, so that we 
can ignore an error introduced by the integrator. On the contrary, using the Floquet matrix 
method we can reliably treat weak electric fields but, since the method explicitly uses the 
boundary conditions (fl3l) . an irrational a has to be approximated by a sequence of rational 
numbers, which makes calculations very time consuming. We did our best employing both 
methods but found no qualitative difference with the case of rational a. This is consistent 

43| where the rationality of a is shown to affect neither dynamics 



with results of Ref. 



42 



nor spectral properties of the system in the tight-binding approximation. 

V. CONCLUSIONS 

We analyzed the problem of the Landau-Zener tunneling for quantum particle in a two- 
dimensional periodic potential, subject to (real or synthetic) in-plane electric and normal 
to the plane magnetic fields. It was found that strong electric field induces transitions 
to higher energy states in the presence of magnetic field as well. Moreover, the estimate 
(jlj) for the rate of LZ-tunneling, which was obtained for zero magnetic field, can be also 
used in the case of a finite magnetic field. This result is of fundamental importance for 
understanding the cyclotron-Bloch dynamics of a quantum particle. In particular, it was 
shown in 43|, |44j devoted to dynamical and spectral properties of the system ([I]) in the tight- 
binding approximation that the particle eigenstates are localized for almost all directions of 
the electric field, where the localization length tends to one lattice site if F tends to infinity. 
The present analysis, however, shows that the tight-binding approximation is valid up to 
some finite F while LZ-tunneling across the energy gap can be neglected. Thus the wave- 
function localization in the limit of large F can be questioned. We reserve this problem for 
future studies. 
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